Entanglement dynamics for qubits dissipating into a common environment 
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We provide an analytical investigation of the entanglement dynamics for a system composed of an arbitrary 
number of qubits dissipating into a common environment. Specifically we consider initial states whose evolution 
remains confined on low dimensional subspaces of the operators space. We then find for which pairs of qubits 
entanglement can be generated and can persist at steady state. Finally, we determine the stationary distribution 
of entanglement as well as its scaling versus the total number of qubits in the system. 
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I. INTRODUCTION 

Entanglement is synonymous of quantum correlations that 
cannot be explained by any local classical theory. Initially 
this notion was relegated to foundational issues, but in the last 
decades it has pervaded quantum information and other areas 
of physics . Being a purely quantum phenomenon it was 
considered fragile under contamination of environment noise 
in open quantum systems. Nevertheless, in recent years it has 
been shown that it can persist for long time up to stationary 
conditions in several contexts J2l (see also ^ for a striking 
experiment on long living entanglement). More importantly it 
has been shown that the environment can play a constructive 
role in establishing entanglement ||4j, |5[] . Quite remarkably it 
happens that even without any interaction among subsystems 
a common dissipative environment is able to induce entangle- 
ment @]. However, this possibility has been proved true only 
for systems composed by two subsystems due to technical dif- 
ficulties arising when accounting for more subsystems. 

Here we overcome these difficulties and achieve results for 
an arbitrary number of subsystems by considering a limited 
number of initially available excitations and by exploiting a 
power expansion for the dissipation superoperator. This al- 
lows us to restrict the analysis to a low dimensional subspace 
of the operators space. 

We find that pairwise entanglement is not created at any 
time between initially excited qubits. Instead, it is created and 
persists at steady state for pairs of initially excited and initially 
not excited qubits. The same holds true for pairs of qubits 
initially in the ground state, but in such a case the amount of 
entanglement is negligibly smaller than the previous case. 



II. SYSTEM DYNAMICS 

Let us consider a system of n qubits with associated Hilbert 
space V. ~ C 2(Sn . Let {|0}, be the orthonormal basis, 

with |0) (resp. |1)) the ground (resp. excited) single qubit 



state. 



We want to study the dynamics of this system of qubits 
when they dissipate together into an environment at zero tem- 
perature. It will be governed by the Lindblad master equation 



p(t) = 2ap(t)o- t - a^ap(t) - p(t)^ a = T>p(t), (1) 



where 



(2) 



with Gi := |0)(1| for the ith qubit. The dissipation rate has 
been set equal to 1 for the sake of simplicity. 

Let us describe the strategy we will put into practice to 



solve fl}. The fomal solution reads p(t) 



a tT> 



p(0). To ex- 



plicitly find p(t) let us first write the Taylor expansion: 



p(t) = p(0) + Wp(0) + yV 2 p(0) + -V 3 p(0) + 



3! 



(3) 



Now notice that once the initial state p(0) has been chosen, 
the dynamics (being purely dissipative at zero temperature) 
can only decrease the number of initial excitation present in 
p(0). That is, repeated applications of V to p(Q) will leave 
the state within a subspace H p ( ) C H of the Hilbert space 
H = n <g> H* where H* stands for the dual of H. After hav- 
ing identified H p ( ), i.e. a set of operators on Ti. spanning 
Hp(o), we will write down p(0) as linear combination of such 
operators with unknown time dependent coefficients. Then, 
by inserting this expansion into Eq.(Q]) we will derive a set of 
linear differential equations for the unknown coefficients. In 
this way if the initial state contains a small number of exci- 
tation we can hope to provide an analytical solution. In fact 
given a number of initial excitations e (e <C n) the following 
inequality holds true: 



dirrM p{0) < (2 e ) 2 < (2") 2 = dimU. 



(4) 
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Finally, notice that ([TJi does not satisfy the condition for the 
uniqueness of a stationary solution. This condition requires 
that the only operators commuting with the Lindblad operator 
a must be multiples of identity 171] . Hence we should expect 
different steady states depending on the choice of p(0). 
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III. ONE INITIAL EXCITATION (e = 1) 

In this section we assume that initially only the fcth qubit is 
in the excited state while all the others are in the ground state. 
Let us denote this state as fc), then p(0) = \k)(k\. 

By applying V to p(0) = \k)(k\ we can easily find the 
following closed relations: 



V\k){k\ 
V\G){G\ 
V(\E ¥ )(k\ + \k)(E ¥ \) 



2\G)(G\-(\E ¥ )(k\ + \k)(E ¥ \), 
0, 

i(n-l)\G){G\-2\E ¥ )(E ¥ \ 
-2(n-l)|fc)(fc| 
-n(\E ¥ )(k\ + \k)(E ¥ \), 
V\E ¥ ){E ¥ \ = 2(n-lf\G){G\ 

-2(n-l)\E ¥ )(E ¥ \ 
-(n-l)(\E ¥ )(k\ + \k)(E ¥ \), 

(5) 



where we have defined 



\G) := |0)« 



\Ey) 



EN) 



(6) 



stands for a state in which the ith qubit is in the excited 
state and all the others are in the ground state). 
With the help of Eq.© we identify the subspace 



I p(0) = S pan{\G)(G\,\k)(k\,\E ¥ )(E ¥ \, 
(\E ¥ )(k\ + \k)(E ¥ \)}. 



(7) 



Then we can write 



p(t) = a {t)\G){G\+ ai {t)\k){k\+a 2 {t)\E ¥ ){E ¥ \ 

+a 3 (t)(\E ¥ ){k\ + \k){E ¥ \), (8) 

which, upon insertion into (Q~|i, leads to the set of differential 
equations 

d = 2a! + 2(n — l) 2 a 2 + 4(n — l)a 3 , 
di = —2a! — 2(n — l)a 3 , 
d 2 = — 2(n — l)a 2 — 2a 3 , 

d 3 = -ai - (n - l)a 2 - na 3l (9) 

with initial conditions ao(0) = a 2 (0) = a 3 (0) = 0, ai(0) = 
1. 

Solving them we arrive at the density operator 

p(t) = (l-f(t)f\k)(k\ + f(t)(2-nf(t))\G)(G\ 
+f(tf\E ¥ )(E ¥ \ 

-f(t)(l-f(t))(\E ¥ )(k\ + \k)(E ¥ \), (10) 



where /(t) :=i(l 



')• 



To analyze the pairwise entanglement between the qubit k 
and a generic other qubit j we compute the partial trace of 




FIG. 1: Concurrence Ck,j versus time t. From top to bottom the 
system size is n = 2, 6, 10, 14. 
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FIG. 2: Stationary concurrence Ck,j versus system size n. 

( fTUt overall the other qubits and obtain the following density 
operator 

p kd = (l-/(t)) 2 |10)(10| + / 2 (t)|01)(01| 

-/(*)(! - /(*)) (|10)<01| + |0L><L0| - 2|00)(00|) . 

(11) 

Using the concurrence |@] as a measure of entanglement of 
Pk.j in Eq.dTTT) we find its time evolution as 



Ck,j(t) 



(12) 



As it is shown in FigQ] entanglement is generated between the 
fcth qubit (initially in the excited state) and an arbitrary qubit j 
(initially in the ground state), although there is no direct inter- 
action between them. The maximum amount of entanglement 
decreases with the system size n, but the concurrence achieves 
its maximum value faster when the system size is larger. 

Letting t going to infinity in Eq.(fT2l we can find the behav- 
ior of stationary entanglement versus the system size n as 



(13) 



This is shown in FigfJ] 

To study the entanglement between any two qubits initially 
in the ground state, we compute the reduced density matrix of 
jth and mth qubits (j, m ^ fc). By referring to Eq.dTQt it is 
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FIG. 3: Concurrence Cj, m versus time t. From top to bottom the 
system size is n = 2, 6, 10, 14. 




FIG. 4: Pictorial representation of the leading stationary correlations 
(edges) among qubits (circles). Red (reps, blue) circles represents 
qubits initially in the excited (resp. ground) state. 



easy to show that 

Pj , m = (l-2/ 2 (i)) |00)(00| 

+/ 2 (t)(|10) + |01))(((10| + (01|). (14) 

The concurrence ^ of pj jm in Eq.dT4b results 

/ 1 - e~ nt \ 2 
C,.(*) = 2(— — ) • (15) 

Its behavior is shown in Fig[3] From Eq.([l~5l) it is easy to see 
that at steady state it is Cj. m (oo) = \, which is negligible 
compared to Ckj{oo) for large n. Therefore for large sys- 
tems, we have a star graph (see FigHI as steady state. 



IV. TWO INITIAL EXCITATIONS (e = 2) 

In this section we assume that initially two qubits, say the 
kth and the Zth, are in the excited state, while all the others 
are in the ground state. Let us denote this state as \k, I), then 
p(0) = \k,l)(k,l\. 

By applying V to p(0) = \k, l)(k, l\ we can find the follow- 
ing closed relations: 



V\k,l)(k,l\ = 


2\k- 


\-l)(k + l\ -A\k,l){k,l\ — A, 


V\k + l)(k + l\ = 


8\G)(G\-2n-4\k + l)(k + l\, 


V\G)(G\ = 


0, 




vn = 


2(n 


- 2) (4\G)(G\ - |fc + l)(k + l\) -nil- 4\Eyi)(E V i}, 


V\E W )(E^ = 


(n- 


2) (2{n-2)\G){G\-Q-2\E^(E^), 


OA = 


4(n 


-2)(\k + l){k + l\ - \k,l)(k,l\)-2\B)(B\+4(Sl-IL) - (n + 2)A, 


V\B)(B\ = 


2(n 


- 2) 2 |fc + Z)(fc + !| + 2(?i - 2)(2fi - A) + 8|%}(£W - 2n\B)(B\ - 


vu = 


(n- 


3) (2f2 - A) - 2(n - 2)n - T, 


VT = 


2(n 


- 3) (ME^iE^ - \B)(B\) + 2(n - 2)(n - 3)fl - 8|fl)<F| - 2(n - 


V\H){H\ = 


(n- 


3) (2(n-3)\E^{E^-4\H)(H\-T), 



4r, 



(16) 



where, in addition to Eq.©, we have defined 



k + l) 


■= \k) + \l), 




\E W ) 


n 

= E ^ 










\B) 


■= El*' fc > 4 










\H) 







h,i ^ k, /, 



h>i 



(17) 



and 

Q := |^(fc + /| + |fc + /)(£^, 
A := |£?)<M| + |M><£|, 
II := |if}(M| + \k,l)(H\, 
T := |B)(JT| + |JT)(B|. (18) 
With the help of Eq.JToTl we can identify the subspace 

H p(0) = span{\G){G\,\E ¥ ) ) {E ¥l ) [ ,\k + l){k + l\,n, 
\k,l}(k,l\,\B}(B\,\H){H\,A,n,T}. 

(19) 
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Then we can write 

P(t) 



following set of differential equations: 



b (t)\G)(G\ + b^E^iE^ + b 2 (t)\k + l)(k + l\ 
+b 3 (t)n + b i (t)\k,l)(k,l\ + b 5 (t)\B)(B\ 

+b 6 {t)\H)(H\ + b 7 (t)A + b B {t)n + b 9 (t)r. (20) 



b (t) = 2{n-2) 2 b 1 {t) 
v(f) = Mv(i), 



i2(t) + 8(n-2)6 3 (*),(21) 
(22) 



By inserting this expansion of p(t) into Eq.(Q]i we arrive at the where v = (pi, b 2l &3 ; &4, £>5, be, 67, 6s, bg) T and 
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(23) 



with the initial conditions &o(0) = 0, v(0) = 
(0,0,0, 1,0, 0,0,0, 0) T . 

It is worth remarking that Eq.(t2Tb has been kept out of 
Eq.(l22l to make the matrix M non-singular. Still M is not 
normal, hence cannot be diagonalized. Nevertheless Eq.(l22b 
can be solved using the Laplace transform 



v(s) 



*v(rj)di 



(s/ 9x9 - M) _1 v(0), 



(24) 



where ^9x9 denotes the 9x9 identity matrix. By inverting the 
Laplace transform we get the coefficients bi, i = 1, ...9. 
Then, bo straightforwardly follows by integrating Eq.(l2ll. 
Their explicit expressions are given in Appendix A. 

The reduced density operator of qubits k and I (initially in 
excited states) results from Eq.(|20b as 



p M (t) = 6 4 (*)|ll)(ll| 

+ (b 2 (t) + (77 - 2)6. (f)) (|10) + |01))((10| + (01|) 



[b (t) + (n-2)[b 1 (t) + ^—^-b 6 (t) 



|00)(00|. 
(25) 



By computing the concurrence ^ of Pk,i(t) it is possible to 
see that the two qubit are not entangled at any time. 

The reduced density operator of qubit k and a generic qubit 



j initially in the ground state derived from Eq.(T20t is given by 

AuW = & 5 (t)|Ll><ll| 

+ (b 2 (t) + (n-3)b 5 (t) + K(t))\W)(W\ 

+ (6 1 (t)+6.(t) + (n-3)6 8 (t))|01)<01| 

+ (6,(t) + (n-3)6 B (t)(|10)<01| + |01>(10|) 

+ (b (t) + (77 - 3)6, (t) + b 2 (t) + (77 - 3)6 5 (t) 

(77 - 3)(77 - 4) 



+ - 



Then, the concurrence 18] of the state 



K(t) |oo)(oo|. 



results 



(26) 



C kJ (t) = -2(b 3 (t) + (n-3)b 9 (t)+b 7 (t)) 

b 5 (t)(b (t) + (n - 3)h(t) + b 2 (t) 

on, , \ (n-3)(n-4), , ,nV2 
+(77 - 3)fe. 5 W + - 1 -& 6 («) • 

(27) 

Blue lines in Fig|5] show the behavior of entanglement be- 
tween qubits k and j quantified by Ck.j vs time. It is generated 
as qubits start interacting with the environment and, as it can 
be seen, the stationary entanglement is achieved faster as the 
system size increases, though the maximum value decreases. 
In the limit of t — > 00 Eq.d27l) becomes 



Ck,j(oo) 



2(n 5 



II77 4 + 3977-* - 53n 2 + 2477 - 



77 2 (77 - 1) 2 (77 - 2) 5 



2(77-3) 



l) 2 (77-2) 2 
1/2 



777 4 - 5077 3 + 119n 2 



10477 + 22 



(28) 
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FIG. 5: Blue lines: concurrence Ck.j versus time t. Green Lines: 
concurrence Cj >m versus time t. From top to bottom the system size 
is n = 4, 7, 10, 13. 




FIG. 6: Pictorial representation of the leading stationary correlations 
(edges) among qubits (circles). Red (reps, blue) circles represents 
qubits initially in the excited (resp. ground) state. 



which for large n behaves like — . 

The other possible case study is the entanglement between 
qubits initially in ground sate, say qubits j and m (j, m ^ 
k, I). Their reduced density operator, derived from Eq.(l20b. 
reads 

p Jim = 6 e (t)|ll)(ll| 

+(&,(*) + 26, (t) + (n-4)6„(t)) 

x(|10) + |01))((10| + (01|) 

+ {b (t) + (n-A)b 1 (t) + 2b 2 (t)+b 4 (t) 



+2(n - 4)6, (t) + ( " 4) , (n 5) b e (t) ) |00) (00| . 



(29) 



The concurrence (s| of the state d29l reads 

C j>m (t) = 2(b 1 (t)+2b 5 (t) + (n-4)b 6 (t)) 

b 6 (t)(b (t) + (n - 4)6i(i) + 2b 2 (t) + b A (t) 

^, , x (ra-4)(n-5), . ^i 1 ^ 
+2(n-4)& 5 (t) + ^ 'i %{t)) . 



(30) 



In the stationary condition C» m becomes 

4(n 4 - 2n 3 - In 2 + lOn - 4) 



(7j, m (oo) 



n 2 (n- l) 2 (n-2) 2 
4 



n(n-l) 2 (n-2) 5 



104n J + 180n z - 152n - 



n 6 - 10n 5 + 41n 4 

1/2 



(31) 



which vanishes for large system size. Green lines in FigO 
show that the amount of entanglement generated between 
qubits initially in the ground state is negligible compared to 
that of qubits k and j. Therefore in the stationary state we 
have a bipartite graph as depicted in Fig|6] 



V. CONCLUDING REMARKS 



In conclusion, we have studied the entanglement dynamics 
of a system composed of an arbitrary number of qubits dissi- 
pating into a common environment. 

The consideration of models like this is becoming increas- 
ingly pressing with the continuing miniaturization of infor- 
mation processing devices. For instance, in solid state im- 
plementations qubits may be so closely spaced that they can 
experience the same environment |@]. 

By considering initial states whose evolution is confined on 
low dimensional subspaces of the operators space we found 
that pairwise entanglement is not created at any time between 
initially excited qubits. Instead, it is created and persists at 
steady state for pairs of initially excited and initially not ex- 
cited qubits. The same holds true for pairs of qubits initially 
in the ground state, but in such a case the amount of entangle- 
ment is negligibly smaller than the previous case. This lead 
to the emergence of structures (subgraphs) which makes this 
subject appealing also from the perspective of complex net- 
works iflOl . 

Finally, it is worth noticing that the extension of the pre- 
sented analysis to states with an arbitrary number e of initial 
excitations can be done in a computationally efficient way due 
to the polynomial scaling of dimM p ^ vs e. We can arrive at 
such a scaling (much tighter with respect to the bound given 
in Eq.©), by simply observing the results in Sec. Ill and IV. In 
fact, in Sec. Ill the density matrix can be written, in the vector 
basis {\G), \k), \E V )}, as 



a 
a\ 03 
a 3 a 2 



(32) 



Furthermore, in Sec. IV the density matrix can be written, in 



6 



the vector basis {\G), \E ¥J ), \k + I), \k, I), \B), \H)}, as 



(33) 



/6 \ 

bi b 3 

6 3 62 

64 b 7 b 8 

b 7 b 5 b 9 

V b 8 b 9 b 6 J 



Then, given the diagonal block structure of the above matri- 
ces (|32| | and d33l . we can argue that for a number of initially 



excitations e the density matrix will have e + 1 blocks of di- 
mensions 1, 2, • • ■ ,e + 1. Each block of dimension r (with 
1 < r < e + 1) has l "( r + 1 ) number of real parameters. Since 
the total number of parameters in the matrix determines the di- 
mension of Hp( ), we will have dimM p ^ = Ylr=\ 1 ^ ? " 2 I ~ 1 ' ) = 
(1 + e)(2 + e)(3 + e)/6 which is polynomial of degree 3 of e. 



Appendix A: Solution of Eqs. {21) and @2) 



-4(n-l)t _ ^ g -2nt _ j 



&0 



6l = _±_(l_ e -™t)2 ( l_ e -2(n-2) t)) 



&2 
&3 
&4 
&5 



n(n-l)(n-2) 2 \(n - l)(n - 2) J ' 

/ e -2(n-l)i e -(ra-2)i 2 

K = 41 — r 2— + - 

b \n(n-l) n(n-2) (n-l)(n-2) 

e -4(n-l)t (g n _ 8 ) e -(3»-4) ( n _ 3 )2 e -2(n-l)t („ _ 4)e -2(n-2)t 

7 ~'n 2 (n-l) 2 + n 2 (n-l)(ri-2) + n(n - l) 2 (n - 2) + ?7 2 (n-2) 
(n-3)(n~6)e-("- 2 )* (n - 3) 2 

+ n(n- l)(n-2) (n - l) 2 (n-2)' 

e _4(„-l)t / _ 4 N e -(3n-4)t 



n - 2 V 4(n - 1) 2n 
4 

(n ^2)n : 
(l-e- 2 ("- 2 )«)(2e-" t + n-2) 2 
(n - 2)n 2 ' 
2(1 - e-" f )(l - e- 2 ("- 2 ) f )(n - 2(1 - e~ nt )) 
(n - 2)n 2 : 

2(n - l)e-("- 2 )* + 2e- 2 ("~ 1 > t + n(n - 3) 
ri(n — 1) 

4e -4(n-l)t 4 ( n 4 ) e -(3n-4)t ( n - 4 )2 g -2(n-2)t 4 ( n _ 3 )g-2(n-l)t 

(n - l) 2 ?i 2 + n 2 (n-l)(n-2) + 7i 2 (n - 2) 2 n(n - l) 2 (n-2) 

2(n-3)(n-4)e-("- 2 )* / (n - 3) 



7i 2 (n-l) 2 n 2 (n- l)(n-2) 

(n 2 -5n + 8)e- 2 ("- 1 ) t e" 2 ^" 2 )* (n - 4) e -(™- 2 )* , (n - 3) 



2n(n- l) 2 (ri- 2) n 2 (n-2) n(n-l)(n-2) 2(n - l) 2 (ri - 2) 

2 ^ 2 n 2 (7i-l) 2 ~ ? i 2 (n-l)(n-2) 

(n-5)e- 2 ("- 1 )* (n- \) e -^-2)t ^ _ 10 ) e -(n-2)t ( ? i - 3) 

V n{n- l) 2 (n-2) + n 2 (n - 2) 2 n(ra - l)(n - 2) 2 + (n - l) 2 (n - 2) 2 1 

I 
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